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STATEMENT  OF  THE  PROBLEM 


Although  most  of  the  work  on  helicopter  rotor  blades  presented  in  the 
literature  to  date  is  based  on  a  set  of  differential  equations  of  motion  that 
include  only  quadratic  nonlinearities  (obtained  by  expanding  the  full 
nonlinear  equations  about  the  undeformed  state  of  the  blade),  it  has  been 
recently  determined  that  cubic  nonlinearities  can  play  a  significant  role  in 
the  equilibrium  and  response  of  a  rotor  blade  in  hover.  This  is  to  be  expected 
since  equations  where  only  quadratic  nonlinearities  are  retained  do  not 
reduce  to  the  proper  nonlinear  equations  of  motion  when  the  blade  is  not 
rotating. 

One  of  the  objectives  of  this  work  was  to  provide  a  thorough  understanding 
of  the  importance  and  treatment  of  nonlinearities  in  the  response  and 
stability  of  helicopter  rotor  blades,  and  to  pinpoint  which  approximations 
caused  by  "ordering  schemes"  lead  to  troblesome  results.  Another  objective 
was  tc  refine  a  perturbation  methodology,  valid  for  the  cases  of  hover  and 
forward  flight,  that  was  initiated  by  the  principal  investigator.  The  analysis 
of  the  response  of  rotor  blades  by  such  a  perturbation  methodology  provides 
an  alternative  to  the  numerically-intensive  Floquet  theory.  It  also  provides, 
with  the  aid  of  symbolic  computation,  closed  form  approximations  for  the 
various  levels  of  approximation  of  the  response,  and  discloses  a  number  of 
important  features  of  the  response  for  different  flight  conditions.  In 
addition,  the  effect  of  truncating  higher  order  nonlinearities  on  the 
equilibrium  state  and  eigenvalues  associated  with  perturbed  motions  of 
beams  and  rotor  blades  was  assessed  in  detail. 
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SUMMARY  OF  MOST  IMPORTANT  RESULTS 


It  is  now  well  known  that  the  nonlinear  differential  equations  of  an  initially 
straight  non-rotating  beam  have  cubic  nonlinearities  that  are  due  to  both 
curvature  and  inertia  effects.  These  nonlinearities  can  play  a  significant  role 
in  the  response  of  the  beam  to  external  loads.  Although  most  of  the  work  on 

helicopter  rotor  blades  presented  in  the  literature  to  date  is  based  on  a  set 

of  differential  equations  of  motion  that  include  only  quadratic 

nonlinearities,  it  has  been  recently  determined  that  cubic  nonlinearities  can 
play  a  significant  role  in  the  equilibrium  and  response  of  a  rotor  blade  in 
hover.  This  is  to  be  expected  since  equations  where  only  quadratic 

nonlinearities  are  retained  do  not  reduce  to  the  proper  nonlinear  equations  of 
motion  when  the  blade  is  not  rotating. 

There  are  a  number  of  analyses  of  rotor  blade  dynamics  in  the  helicopter 
industry.  The  computer  programs  for  those  analyses  predict  different 
eigenvalues  associated  with  the  perturbed  motion  of  the  blade.  Also,  in  some 
of  the  published  work  available  in  the  open  literature,  there  have  been  a 
number  of  questions  regarding  the  proper  treatment  and  influence  of 
nonlinear  terms,  especially  quadratic  and  cubic  order  terms,  in  the  dynamics 
of  a  rotor  blade.  These  problems  were  addressed  in  this  research  work  and 
the  effect  of  truncating  higher  order  nonlinearities  on  the  equilibrium  state 
and  eigenvalues  associated  with  perturbed  motions  of  a  non-rotating  beam 
and  of  a  rotor  blade  were  assessed  in  detail. 

in  our  work,  special  care  was  taken  in  formulating  the  differential  equations 
of  motion  for  a  blade  in  order  to  be  able  to  identify  which  commonly  used 
approximations  lead  to  invalid  results.  We  have  been  able  to  show  [1]  that 
approximating  all  the  diagonal  elements  of  transformation  matrices  by  unity 
while  retaining  nonlinear  terms  off  the  diagonal  is  one  of  such  erroneous 
approximations.  This  approximation  alone  renders  the  subsequent  nonlinear 
equations  invalid  for  analyzing  the  dynamics  of  beams  or  rotor  blades. 
Equations  obtained  when  such  approximation  is  used  lead  to  results  that 
depend  on  Euler  angle  sequence,  which  is  an  incorrect  result  even  on  physical 
grounds  alone.  In  developing  this  work  the  analysis  of  a  cantilever  beam  with 
a  tip  mass  proved  to  be  invaluable.  Complete  details  of  this  work  were 
reported  in  [1].  The  work  was  done  using  the  facilities  of  our  Symbolic 
Computation  Laboratory.  We  made  use  of  MACSYMA  to  develop  and  to  analyze 
the  equations  of  motion  and  to  generate  error  free  extensive  Fortran  code 
which  was  needed  to  generate  all  the  graphical  results  presented  in  [1],  The 
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effect  of  analyzing  the  dynamics  of  a  beam  by  replacing  the  full  nonlinear 
differential  equations  of  motion  by  a  set  of  nonlinear  equations  expanded 
about  the  undeformed  state  of  the  beam,  when  such  state  is  not  an 

equilibrium,  was  studied  in  detail  [2],  In  addition,  the  use  of  a  transfer 
matrix  methodology,  instead  of  Galerkin's  method,  was  also  explored  in  [2], 
Based  on  the  results  and  conclusions  presented  in  [2],  the  differential 

equations  of  motion  for  a  rotor  blade  were  reformulated  and  the  blade's 
response  in  hover  was  analyzed  with  such  methodology.  For  this,  an 
essentially  exact  equilibrium  solution  for  the  blade  was  obtained  by 
numerically  integrating  its  full  nonlinear  differential  equations  of  motion. 
For  the  case  of  hover,  considered  in  [3],  the  differential  equations  governing 
the  equilibrium  state  of  the  blade  are  a  set  of  nonlinear  ordinary  differential 
equations.  Standard  IMSL  routines  were  used  to  solve  the  two-point  boundary 
value  problem  represented  by  such  equations.  When  the  equilibrium  state  is 
perturbed,  and  the  resulting  equations  are  linearized  for  infinitesimally 
small  perturbations,  the  solution  to  the  linearized  equations  can  be  written 
as  E(s)  erf  where  s  denotes  distance  measured  along  the  blade  and  t  denotes 

time;  JF(s)  is  a  column  of  state  variables  that  were  properly  chosen  so  that 

the  perturbed  differential  equations  are  written  as  a  set  of  first  order 
equations  that  can  be  integrated  numerically  [3].  The  resulting  differential 
equations  for  F(s)  contain  the  complex  quantity  r  as  a  parameter.  Those 
equations  were  integrated  in  the  same  manner  described  in  [2]  to  determine 
the  eigenvalues  r.  Plots  showing  how  the  equilibrium  solution  and  the  real 
and  imaginary  parts  of  r  depend  on  a  number  of  important  parameters  for  the 
system  were  presented  in  [3].  The  results  obtained  were  compared  with  those 
generated  by  using  a  Galerkin  procedure  with  non-rotating  beam  modes  and 
significant  differences  were  noted  for  several  cases.  Such  results  could 
serve  as  a  gauge  for  validating  the  approximate  results  obtained  with 
Galerkin's  procedure. 

It  should  be  noted  that  by  avoiding  Galerkin's  procedure,  such  as  done  in  the 
work  presented  in  [2,3],  the  calculation  of  a  large  number  of  Galerkin 
coefficients  was  avoided  and  the  use  of  ordering  schemes  was  not  necessary. 
Also,  the  eigenvalues  associated  with  the  perturbed  motion  were  calculated 
without  the  need  for  iterating  on  the  number  of  trial  modes,  such  as  non¬ 
rotating  beam  modes,  required  by  the  Galerkin  procedure.  Another  advantage 
of  the  me  work  presented  in  [2,3]  is  that  the  solution  obtained  for  the 
equilibrium  state  the  rotor  blade  was.  waco^i  for  small  :cund-cff  enuis 
inherent  in  any  numerical  scheme,  essentially  exact  since  the  full  nonlinear 
differential  equations  were  used  in  the  numerical  integration.  Such  solution 
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could  also  have  been  obtained  by  other  numerical  methods,  such  as  finite 
elements,  without  major  difficulties.  Additional  details  of  part  of  the  work 
reported  in  [1]  are  included  in  [4,5]. 

Another  objective  of  this  work  was  to  refine  the  development  of  a 
perturbation  methodology  valid  for  the  cases  of  hover  and  forward  flight.  In 
contrast  with  Floquet  theory,  a  perturbation  analysis  provides  significant 
insight  into  the  physics  of  the  problem.  In  a  perturbation  analysis  one  wants 
to  render  the  first  order  differential  equations  linear  and  to  make  the 
analysis  valid  for  any  flight  condition. 

The  perturbation  methodology  is  based  on  the  method  of  multiple  time 
scales,  and  is  valid  both  for  the  cases  of  hover  and  of  forward  flight.  In  order 
to  make  the  perturbation  analysis  valid  for  both  of  these  cases,  the  advance 
ratio  cannot  be  used  as  the  perturbation  parameter.  Instead,  a  small  arbitrary 
"bookkeeping  parameter"  that  is  introduced  artificially,  and  that  linearizes 
the  differential  equations  of  motion,  was  used.  The  significant  amount  of 
"algebra"  associated  with  the  perturbation  analysis  was  done  with 
computerized  symbolic  manipulation  using  MACSYMA. 

It  should  be  noted  that  perturbation  analyses  have  been  used  by  a  few  other 
researchers  in  helicopter  dynamics.  However,  the  analyses  that  have  been 
presented  in  the  literature  to  date  deal  separately  with  the  cases  of  hover 
and  of  forward  flight.  For  each  case  a  different  small  perturbation  parameter 
was  used,  such  as  the  collective  pitch  angle  for  the  case  of  hover  and  the 
advance  ratio  for  the  case  of  forward  flight.  However,  none  of  these 
parameters  render  the  equations  of  motion  linear  in  the  limit  as  they 
approach  zero.  Thus  such  parameters  are  not  ideally  suited  for  a  perturbation 
analysis  in  the  presence  of  nonlinearities. 

In  order  to  make  our  perturbation  analysis  valid  for  both  cases  of  hover  and 
of  forward  flight,  and  also  to  make  it  handle  nonlinearities,  a  small  arbitrary 
"bookkeeping  parameter"  is  introduced  in  the  analysis.  As  in  other  research 
work  published  by  the  principal  investigator,  such  parameter  is  used  to  keep 
track  of  orders  of  smallness  from  a  perturbation  point  of  view.  It  should  be 
noted  that  ordering  here  does  not  have  the  same  meaning  as  that  of  ordering 
schemes  that  have  been  used  elsewhere  to  eliminate  complicated  terms  in 
the  analysis. 

The  details  of  the  perturbation  analysis,  including  the  effects  of  nonlinear 
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terms  up  to  0(e3),  are  described  in  [6].  An  alternative  perturbation  analysis 
based  on  a  set  of  differential  equations  used  for  a  Floquet  stability  analysis 
was  also  performed.  The  latter  analysis  is  simpler  than  the  full  nonlinear 
perturbation  analysis.  It  was  found  that  the  two  analyses  produce  essentially 
the  same  results  for  advance  ratios  up  to  about  0.3,  which  is  in  the  range  of 
practical  interest  for  most  applications.  The  perturbation  analysis  also 
dicrioses  the  resonance  conditions  in  the  system,  clearly  points  to  them  as 
the  cause  for  instabilities,  and  allows  one  to  concentrate  the  search  for 
unstable  regions  in  the  system's  parameter  space  near  those  lines  defined  by 
the  resonance  conditions. 

The  refined  perturbation  methodology  developed  in  this  work  consists  in 
first  determining  a  particular  solution  to  the  full  nonlinear  differential 
equations  of  motion.  This  was  done  numerically.  The  particular  solution  is 
either  an  equilibrium  solution,  for  the  case  of  hover,  or  a  periodic  solution 
that  satisfies  some  imposed  conditions.  The  latter  is  the  "trim  solution"  for 
the  case  of  forward  flight.  By  transferring  the  periodic  coefficient  and  the 
damping  terms  from  the  0(e)  to  the  next  order  approximation,  solvability 
conditions  (i.e.,  conditions  for  elimination  of  secular  terms)  were  obtained  at 
both  the  0(e2)  and  at  the  0(e3)  approximations.  To  obtain  the  solvability 
conditions  one  has  to  distinguish  between  several  resonant  conditions 
disclosed  by  the  perturbation  analysis  for  the  0(e2)  and  0(e3)  differential 
equations.  For  the  case  of  a  hinged  spring-restrained  rigid  blade,  it  was 
found  that  the  perturbation  analysis  discloses  the  existence  of  two 
important  resonance  conditions,  namely  0)^=  0^  and  a>^=  1  +  to^.  Here,  co^  and  co^ 
denote,  respectively,  the  lead-lag  and  the  flap  natural  frequencies.  The 
solvability  conditions  for  the  motion  "near"  these  resonance  conditions  were 
obtained  and,  from  these,  the  stability  of  the  motion  was  then  determined. 
Since  the  0(e")  ,  n>1,  differential  equations  obtained  from  the  perturbation 
analysis  are  equations  with  constant  coefficients,  stability  determination 
using  such  equations  is  readily  done  using  the  Routh-Hurwtiz  criterion. 
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